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After a short review of the ultrahigh energy cosmic ray puzzle - the apparent obser- 
vation of cosmic rays originating from cosmological distances with energies above the 
expected Greisen-Zatsepin-Kuzmin cutoff 4 X 10 19 eV - we consider strongly interacting 
neutrino scenarios as an especially interesting solution. We show that all features of 
the ultrahigh energy cosmic ray spectrum from 10 17 eV to 10 21 eV can be described to 
originate from a simple power-like injection spectrum of protons, under the assumption 
that the neutrino-nucleon cross-section is significantly enhanced at center of mass ener- 
gies above fa 100 TeV. In such a scenario, the cosmogenic neutrinos produced during the 
propagation of protons through the cosmic microwave background initiate air showers in 
the atmosphere, just as the protons. The total air shower spectrum induced by protons 
and neutrinos shows excellent agreement with the observations. We shortly discuss TeV- 
scale extensions of the Standard Model which may lead to a realization of a strongly 
interacting neutrino scenario. We emphasize, however, that such a scenario may even 
be realized within the standard electroweak model: electroweak instanton/sphaleron 
induced processes may get strong at ultrahigh energies. Possible tests of strongly inter- 
acting neutrino scenarios range from observations at cosmic ray facilities and neutrino 
telescopes to searches at lepton nucleon scattering experiments. 



1. Introduction 

The Earth's atmosphere is continuously bombarded by cosmic particles ("rays"). 
Their measured flux extends over many orders of magnitude in energy (cf. Fig. [T]l. 
At energies above 10 15 eV, they are observed in the form of extensive air show- 
ers (EAS's), initiated by inelastic scattering processes of cosmic particles off at- 
mospheric nucleons. Ground-based observatories have measured EAS's with en- 
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ergies up to £< 3 x 10 20 eV, corresponding to center-of-mass (CM) energies 
y/s = y / 2m p E < 750 TeV, where m p is the proton mass. Therefore, the high- 
est energy cosmic rays probe physics beyond the reach of the (Very^) Large Hadron 
Collided ((V)LHC), with a projected CM energy of 14 (200) TeV. In this context, it 
is interesting that the measured cosmic ray flux at the highest energies, E > 10 20 eV, 
represents a puzzle. What is this puzzle about? 
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Figure 1. Compilation of measurements of the flux of cosmic rays. The dotted line shows anE 3 
power-law for comparison. Approximate integral fluxes (per steradian) are also shown (adapte 
from Ref0). 



It hinges on the circumstantial evidence that the cosmic rays above 10 17 5 : 18 5 eV 
originate from cosmological distances (for a recent review, see RefEl). This evidence 
is largely based on the apparent large-scale isotropy in the arrival directions of 
cosmic rays (cf. Fig.EJ). Moreover, whereas there are only very few - if any - nearby 
source candidates, plausible astrophysical sources are most likely to be found only 
at cosmological distances. 

If the highest energy cosmic rays are nucleons (or nuclei), if their sources are 
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Figure 2. Arrival directions of cosmic rays detected by the AGASA and Akeno (A20) experiments 
in equatorial coordinates. Open circles and open squares represent cosmic rays with energies 
(4 -r 10) X 10 19 eV, and > 10 20 eV, respectively. The galactic and super-galactic planes are shown 
by the red and blue curves, respectively. Large shaded circles indicate event clusters within 2.5°. 
The shaded regions indicate the celestial regions excluded by a zenith angle cut of < 45 °. UpdatP 
(June 24, 2003) of the published data from RefH 



indeed uniformly distributed at cosmological distances, and if their injection spec- 
tra are power-laws in energy - a reasonable assumption, in view of the measured 
spectrum in Fig. ^ which appears to be approximately of (broken) power-law type 
over many order of magnitude in energy - then their total flux arriving at Earth 
should show a pronounced drop above the Greisen-Zatsepin-KuzmirP (GZK) "cut- 
off" Eqzk = 4 x fO 19 eV. This is due to the fact that, above this energy, the 
universe becomes opaque to high energy nucleons (and nuclei), due to inelastic 
hadronic scattering processes with the cosmic microwave background (CMB) pho- 
tons. The GZK cutoff is, however, not seen in the data, at least not in a significant 
manner (cf. Fig. EJ- Correspondingly, the events above 10 20 eV in Fig. should 
originate from small distances below 50 Mpc, the typical interaction length of nu- 
cleons above Eqzk- However, no source within a distance of 50 Mpc is known in 
the arrival directions of the post-GZK events a . The basic puzzle is: if the sources of 
ultrahigh energy cosmic rays are indeed at cosmological distances, how could they 
reach us with energies above 10 20 eV? 

At the relevant energies, among the known particles only neutrinos can prop- 
agate without significant energy loss from cosmological distances to us. It is this 
fact which led, on the one hand, to scenarios invoking hypothetical - beyond the 
Standard Model - strong interactions of ultrah igh energy cosmic neutrinoJ^ and, 
on the other hand, to the Z-burst scenark P^ I ^ I 

In the latter, ultrahigh energy cosmic neutrinos (UHECz/s) produce Z-bosons 
through annihilation with the relic neutrino background from the big bang. On 



a The dominant radio galaxy M87 in the Virgo cluster, at a distance of about 20 Mpc, has been a 
source candidate for a long timEU The major difficulty with this idea is the isotropy of the arrival 
distribution. It migh t be overcome by invokin g a. particular galactic magnetic field originating from 
a "galactic wind'>13 Criticisms of this modeBSl have been addressed in 
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Figure 3. Ultrahig h e nergy cosmic ray data with their statist ical errors (tog: combination of 
AkenJHl and AGASy»ED data; bottom: combination of Fly's EytM and HiReJHl data) and the 
predictions arising from a power-law emissivity distribution Q corresponding to sources which 
are uniformly distributed at cosmological distances. The best fits between E— = 10 17 ' 2 eV and 
= 10 20 eV are given by the solid lines and correspond to the indicated values of the parameters 
a and n in the source emissivity distribution. The 2-sigma variations corresponding to the minimal 
(dotted) and maximal (dashed) fluxes are also shown. Other parameters of the analysis were 
£ max = 3 X 10 21 eV, z mi „ = 0.012, and 2 max = 2. From RefE3 



Earth, we observe the air showers initiated by the protons and photons from the 
hadronic de cays of these Z-bosons. Though the re quire d ultrahigh energy cosmic 
neutrino fhrd^is smaller than present upper bounds^!, it is not easy to conceive a 
production mechanism yielding a sufficiently l arge one. In the near fu tur e, UHECV 
detectors, such as the Pierre Auger Observatory^, IceCubJ23l ANITaPH, EUScES 
OWL 26 , and SalSA 27 can directly confirm or exclude this scenari 

Scenarios based on strongly interacting neutrinos, on the other hand, are based 
on the observation that the flux of neutrinos originating from the d ecay of the 
pions produced during the propagation of nucleons through the CMEM^ ^^M^l l 
- the cosmogenic neutrinos - shows a nice agreement with the observed ultrahigh 
energy cosmic ray (UHECR) flux above Eqzk- Assuming a large enough neutrino- 
nucleon cross-section at these high energies, these neutrinos could initiate extensive 
air showers high up in the atmosphere, like hadrons, and explain the existence of 
the post-GZK event This large cross-section is usually ensured by new types of 
TeV-scale interactions beyond the Standard Model, such as arising through gluonic 
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bound state leptond^l through TeV-scale grand unification with leptoquarkJ^, 
through Kaluza-Klein modes from compactified extra dimensionJ^ (see, howe ver, 
Ref !""|), or through p-brane production in models with warped extra dimensions^! 
(see, however, Ref!^); for earlier and further proposals, see RefPI and RefPES 
respectively. 

In this review, we discuss strongly interacting neutrino scenarios as a possible 
solution to the GZK puzzle. We present a detailed statistical analysis of the agree- 
ment between observations and predictions from such scenarios^! Moreover, we 
emphasize an example which - in contrast to previous proposals - is based entirely 
on the Standard Model of particle physics. It exploits non-perturbative electroweak 
instanton-induced processes for the interaction of cosmogenic neutrinos with nuclc- 
ons in the atmosphere, which may have a sizeable cross-section above a threshold 
energy E t h — 0((4nmw / 'aw) 2 ) I '(2m p ) = O(10 18 ) eV, whe re mw d enotes the W- 
boson mass and aw the electroweak fine structure constant! ' ' ^ \ 

Our scenario is based on a standard power-like primary spectrum of protons 
injected from sources at cosmological distances. After propagation through the 
CMB, these protons will have energies below Eqzk, so they can well describe the low 
energy part of the UHECR spectrum. The cosmogenic neutrinos interact with the 
atmosphere and thus give a second component to the UHECR flux, which describes 
the high energy part of the spectrum. The relative normalization of the proton and 
neutrino fluxes is fixed in this scenario, so the low and high energy parts of the 
spectrum are explained simultaneously without any extra normalization. Details of 
this analysis can be found in Refill 

The structure of this review is as follows. In the next section, we review our 
procedure to infer the fluxes of protons and cosmogenic neutrinos at Earth, from an 
assumed injection spectrum at the sources. In Sect.03 various possibilities, including 
the one exploiting electroweak instantons, for a large neutrino- nucleon cross-section 
at high energies are discussed, and the induced air shower rate is calculated. In 
Sect. QJ we present a comparison of the predictions with the observations and a 
determination of the goodness of the fit. Possible further tests are mentioned in 
Sect. while conclusions are given in Sect. El 



2. Proton and cosmogenic neutrino fluxes 

Our analysis^ is based on the assumption of a power-law emissivity distribution 
corresponding to uniformly distributed sources. The emissivity is defined as the 
number of protons per co-moving volume per unit of time and per unit of energy, 
injected into the CMB with energy Ei and characterized by a spectral index a and 
a redshift (z) evolution index n, 

C p = jo Er a (1 + z f 6(E max - E t ) 6(z - z min ) 0(z max - z) . (1) 

Here, jo is a normalization factor, which will be fixed by the observed flux. The 
parameters i? max and z min / max have been introduced to take into account certain 
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possibilities such as the existence of a maximal energy, which can be reached through 
astrophysical accelerating processes in a bottom-up scenario, and the absence of 
nearby/very early sources, respectively. Our predictions are quite insensitive to the 
specific choice for E max , z m ; n , and z maXl within their anticipated values. The main 
sensitivity arises from the spectral parameters a and n, for which we determine the 
1- and 2-sigma confidence regions in Sect. 01 

The propagation of particles can be describe c $U | 44 | 45 | ^ y p b ^ z ^Ei\E) func- 
tions, which give the expected number of particles of type b above the threshold 
energy E if one particle of type a started at a redshift "distance" z with energy 
Ei. With the help of these propagation functions, the differential flux of protons 
(b = p) and cosmogenic neutrinos (b = z^, v{) at Earth can be written as 

Fb{E) -^J Q H(z) dE lTz~' (2) 

In our analysis, we took z max = 2 (cf. Ref.^j), while we choose z m - m = 0.012 in order 
to take into account the fact that within 50 Mpc there are apparently no astrophys- 
ical sources of UHECR's. We used the expression H 2 (z) = Hq [Qm (1 + z) 3 + Q\\ 
for the relation of the Hubble expansion rate at redshift z to the present one. Un- 
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Figure 4. Predicted cosmogenic neutrino fluxes per flavor, F yi +Fy l , £ = e, fi, r, originating from 
a power-law proton source emissivity distribution Q and corresponding to the predicted UHECR 
fluxes in Fig. |3] The "best" predictions for the neutrino spectra are given by the solid lines. The 
2-sigma variations corresponding to the minimal (dotted) and maximal (dashed) fluxes are also 
shown. The dotted band labelled by Auger represents the expected sensitivity of th e Pierre Auge r 
Observatory to u T + v T , corresponding to one event per year per energy decadJ^I From RefffS 
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certainties of the latter, H Q = h 100 km/s/Mpc, wittPZI/i = (0.71 ± 0.07) xj;^, are 
included. CIm and JIa, with Q m + = 1, are the present matter and vacuum en- 
ergy densities in terms of the critical density. As default values we choose £Im = 0.3 
and ft a — 0.7, as favored today. Our results turn out to be rather insensitive to the 
precise values of the cosmological parameters. 

We calculated Pu a (z,Ei;E) in two steps, i) First, the SOPHIA Monte-Carlo 
prograrrPSI was used for the simulation of photohadronic processes of protons with 
the CMB photons. For e + e~ pair production, we used the continuous energy loss 
approximation, since the inelasticity is very small (« 10 -3 ). We calculated the P^ a 
functions for "infinitesimal" steps as a function of the redshift z. ii) We multiplied 
the corresponding infinitesimal probabilities starting at a redshift z down to Earth 
with z = 0. The details of the calculation of the P^ a (z, Ef, E) functions for protons, 
neutrinos, charged leptons, and photons will be published elsewhere^! 

Since the propagation functions are of universal usage, we decided to make the 
latest versions of —dPb\ a /dE available for the public via the World- Wide- Web URL 
www.desy.de/~uhecr . 

As an illustration of the outcome of our propagation codes, we display in Fig. 
the predictions for the proton flux at Earth, originating from a power-like source 
emissivity distribution (JJJ with particular a,n, . . . values indicated on the figure 
and in its caption. A nice fit of the data can be obtained apparently for energies 
below <4x 10 19 eV = Eqzk - more on this in Sect.0J The associated predicted 
cosmogenic neutrino flux, for the same parameter values, is displayed in Fig. 



3. Spectrum of neutrino-induced air showers 

The main assumption of strongly interacting neutrino scenarios is that the neutrino- 
nucleon cross-section er*~^ suddenly becomes much larger than rs 1 mb above center 
of mass energies y/s rj 100 TeV. In this case, the corresponding neutrino interac- 
tion length \ u = m p /al%, with a% = a% + cr™^, falls below X a = 1031 g/cm 2 
- the vertical depth of the atmosphere at sea level - above the neutrino threshold 
energy w 10 19 eV. Here af, c N and cr"^ denote the charged current and the new 
contribution to the cross-section. Above the neutrino threshold energy, the atmo- 
sphere becomes opaque to cosmogenic neutrinos and most of them will end up as 
air showers. Quantitatively, this fact can be described by 

F' V {E) = F U (E) l- e -— = F„(E) x <^ K{E) , (3) 

L J { 1 for \ v {E)<^X{e) 

which gives the spectrum of neutrino-initiated air showers, for an incident cosmo- 
genic neutrino flux F v = J^il^i + ^K] horn Eq. in terms of the atmospheric 
depth X(8), with 9 being the zenith angle. 

Such suddenly increasing cross-sections have been proposed in vario us mod- 
els involving physics beyond the Standard ModeP^I ^^^^^^^^^^^ I ^ . Among 
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Figure 5. Total cross section o(i>N — ► p— brane) in a model with n = 6 extra dimensions, out 
of which m have a size L/L„ = 0.25 in terms of the fundamental Planck length L» = , for 

m = 0, . . . , n— 1 from below. The fundamental Planck scale Afjj = [(2-7r) n /87r] 1 /( n + 2 )M < , has been 
chosen as Mjj = M™ ln = 1 TeV in terms of the minimum p-brane mass Af™ ln . The Standard 
Model charged current cross-section <j(vN — ► £X) is also shown (dotted). From RefE3 



the usual sus pects are TeV-scale gravity scenarios with large or warped extra 
dimensions^! In those, the neutrino-nucleon cross-section may be greatly enhanced 
compared to the Standard Model one. As an exam ple, we demonstrate in Fig.[S]that 
p-brane production in neutrino-nucleon scatterin g I I may reach a cross-section 
of 10 mb at ss 10 19 eV, depending on the parameters of the model. This is in 
contrast to microscopic black hole (= 0-brane) productiorl^l which has generically 
too small a cross-sectiorJ^ to solve the GZK puzzle, within the allowed parameter 
ranges. 

In Fig. |HJ we show another example for a strong enhancement in the neutrino- 
nucleon cross-section, which is based entirely on the Stand ard Mode l, exploit- 
ing non-perturbative electroweak instanton-induced processeJH™ Accord- 
ing to the estimates presented in Fig. H3 the electroweak instanton-induced 
neutrino-nucleon cross-section appears to have a threshold-like behavior at E t h — 
0({Airmw I 'ctw) 2 ) I '(2m p ) — O(10 18 ) eV, above which it quickly rises above 1 mb. 
Our quantitative analysis in Ref.^ was based on the cross-section from Ref@Sl 
(solid line in Fig.EJ, however it is quite insensitive to the exact form of it as long 
as it rises abruptly far above 1 mb. Note that such a behaviour is co nsisten t with 
present upper bounds on electroweak instanton-induced cross-sectionJ^EZI How- 
ever, it is fair to say that there are substantial uncertainties in the predictions in 
Fig. HO the absolute size of the cross-section above the threshold energy may well 
be unobservably small. 
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4. Comparison with UHECR data 

The predicted air shower rate induced by protons and neutrinos is given by 

FprediE; a,n, E mgx , z min , z max , jo) = F P (E; . . .) + F' V (E\ . . .) . (4) 

In RefEHl 

we performed a statistical analysis to compare the prediction , within 
the electroweak instanton scenario from Fig. (solid), with the observations and 
presented a measure for the goodness of the scenario. We gave the best fit to the 
observations and the 1- and 2-sigma confidence regions in the (a,n) plane. 

To start the analysis, we had to convert the measured fluxes, which UHECR 
collaborations usually publish in a binned form, into event numbers in each bin. 
We used the most recent results of the HiRes and AGASA collaborations and did 
our analysis separately with both data sets. We concentrated on the energy range 
IO 17,2 eV - 10 21 eV which is divided into 38 equal logarithmic bins. In the low 
energy region, there are no published results available from AGASA and only low 
statistics results from HiRcs-2. Therefore, we included the results of the predecessor 
collaborations - AkencEl and Fly's Eye, respectively - into the analysis. With a 
small normalization correction, it was possible to continuously connect the AGASA 
dateP with the Akeno ones and the HiRes- 1 monocular data-- with the Fly's Eye 
stereo onei^, respectively (cf. Figs. E3 and EJ). 

The goodness of the scenario was determined by a statistical analysis. We 
determined the compatibility of different (a,n) pairs (cf. Eqs. Q and |@}) with 
the experimental data. For some fixed (a,n) pair, the expected number of events in 
individual bins are (A = {Ai, A r } with r being the total number of bins (in our 
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Figure 6. Predictions of the electroweak instanton-induced neutrino- nucleoli cross-section u ' v 
(solicBHl and dashecEl) in comparison with the charged current cross-section cr cc N (dotted) from 
RefP as a function of the neutrino energy E v in the nucleon's rest frame. 
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Figure 7. Ultrahigh energy cosmic ray data (Akeno + AGASA on the upper panel and Fly's Eye 
+ HiRes on the lower panel) and their best fits (solid) within the electrowcak instanton scenario 
from Fig. [H] (solid) for i? m ax — 3 X 10 22 eV, ^min — 0.012, z m ax — 2, consisting of a proton 
component (dotted) plus a cosmogenic neutrino-initiated component (dashed). From RefESl 



case 38). The probability of getting an experimental outcome k = {k\, ...k r } (where 
ki are non-negative integer numbers) is given by the probability distribution function 
P (k), which is just the product of Poisson distributions for the individual bins. 
We also included the ss 30% overall energy uncertainty into the P(k) probability 
distribution. We denote the experimental result by s = {si, ...s r }, where the Sj-s are 
non-negative, integer numbers. The (a,n) pair is compatible with the experimental 
results if 



£ ^(k)<c. 



(5) 



k|P(k)>P(s) 



For a l-(or 2-)sigma compatibility one takes c=0.68 (or c=0.95), respectively. The 
best fit is found by minimizing the sum on the left hand side. 

Figure0shows our best fits for the AGASA and for the HiRes UHECR data. The 
best fit values are a = 2.68(2.68) and n = 2.65(2.9), for AGASA(HiRes), within the 
electroweak instanton scenario from Fig. [S] (solid). We can see very nice agreement 
with the data within an energy range of nearly four orders of magnitude. The fits 
are insensitive to the value of E max as far as we choose a value above «3x 10 21 eV. 
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Figure 8. Confidence regions in the a— n plane for fits to the Akeno + AGASA data (2-sigma 
(long dashed)) and to the Fly's Eye + HiRes data (1-sigma (solid); 2-sigma (short-dashed)), 
respectively, within the electroweak instanton scenario from Fig. [U] (solid), for E max >3x 10 21 eV, 
Zmin > 0, 2 max = 2. From RefHH 

The shape of the curve between 10 17 eV and 10 19 eV is mainly determined by the 
redshift evolution index n. At these energies the universe is already transparent for 
protons created at z f=a 0, while protons from sources with larger redshift accumulate 
in this region. The more particles are created at large distances - i.e. the larger n is 
- the stronger this accumulation should be. In this context, we note that the data 
seem to confirm our implicit assumption that the extragalactic uniform UHECR 
component begins to dominate over the galactic one already at « 10 17 eV. If we, 
alternatively, start our fit only at 10 18 5 eV - corresponding to the assumption that 
the galactic component dominates up to this energy - we find, however, also a very 
good fit, with a very mild dependence on n and the same best fit values for a, 
with a bit larger uncertainties. The peak around 4 x 10 19 eV in Fig. shows the 
accumulation of particles due to the GZK effect. Neutrinos start to dominate over 
protons at around 10 20 eV. 

It is important to note that, if we omit the neutrino component, then the model 
is ruled out on the 3-sigma level for both experiments. This is due to the fact that we 
excluded nearby sources by setting z m - ul ^ (see also RefP^). The choice z m - ul = 
makes the Hi Res da ta compatible with a proton-only scenario on the 2-sigma level 
(see also RefsP2EHl). 

Figure |H1 displays the confidence regions in the (a,n) plane for AGASA and 
HiRes. The scenario is consistent on the 2-sigma level with both experiments. 
For HiRes, the compatibility is even true on the 1-sigma level. It is important to 
note that both experiments favor the same values for a and n, demonstrating their 
mutual compatibility on the 2-sigma level (see also Ref P"™) . If we ignore the energy 
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uncertainty in the determination of the goodness of the fit, they turn out to be 
inconsistent. 




1 au 1 pc 1 kpc 1 Mpc 

log (size, km) 

Figure 9. The Hillas diagram showing size and magnetic field strengths of possible sites of 
particle accele rati on. Objects below the diagonal lines (from top to bottom), derived from the 
Hillas criteriorl^U E ma x ~ 2Ze B r for the maximum energy acquired by a particle of charge Ze 
traveling in a medium of size r with a magnetic field B, cannot shock accelerate protons above 
10 21 eV, above 10 20 eV and iron nuclei above 10 20 eV, respectively. (This version of the picture 
is courtesy of Murat Boratav). 

Finally, let us emphasize that the same fit results are valid for all strongly in- 
teracting neutrino scenarios, as long as the neutrino-nucleon cross-section has a 
similar threshold- like behavior as in Figs. and with a neutrino threshold energy 
< 4 x 10 19 eV and a cross-section > 1 mb above threshold. It is also important to 
note that the energy requirements on the sources of the primary protons are com- 
paratively mild. To obtain a good fit, we need E max > 3 x 10 21 eV. An inspection 
of the Hillas diagram in Fig. 03 reveals that there are a number of reasonable astro- 
physical source candidates, notably neutron stars and gamma ray bursters (GRB's), 
which may provide the necessary conditions to accelerate protons to the required 
energies by conventional shock acceleration. 
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5. Further tests 

There are a number of further possible tests of strongly interacting neutrino sce- 
narios, ranging from astroparticle tests, which include searches at EAS arrays and 
neutrinos telescopes, to laboratory tests at present and future accelerators. We will 
review some of those in this Section. 



5.1. Astroparticle tests 
5.1.1. Searches at EAS arrays 

One possibility to test the ultrahigh energy neutrino component in the EAS data 
is to study the zenith angle dependence of the events in the lfj 18 ^ 20 eV range, 
which will reflect the energy dependence of the neutrino-nucleon cross-sectiorP^E2. 
Near the threshold energy in strongly interacting neutrino scenarios, there will be 
always a range of energies where the cross-section is already sizable, but does not 
yet reach hadronic values (cf. Figs. |S] and [fJJl, in particular, where c*^ < 0.56 mb, 
corresponding to the atmospheric depth at larger zenith angles, 9 > 70°. Therefore, 
for these energies, neutrino-initiated showers can be searched for at cosmic ray 
facilities by looking for quasi-horizontal air showerd^l > 70°. We have checked 
in RefEOI that the rate from our electroweak instanton prediction (cf. Fig. [(J) is 
consistent with observational constraints found by the Fly's Eye*^ and AGASi^P^l 
collaborations. For the case of p-brane production, such cons traints can be avoided 
in warped extra dimension scenarios with fine-tuned size 
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Figure 10. AGASA integrated energy spectrum of cosmic rays (cl osed circles) and their contri- 
bution to clustered events (open squares, see also Fig.|5J. From RefEU 
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The arrival directions of the cosmogenic neutrinos should pretty much coincide 
with the direction of the primary protons. Therefore, strongly interacting neutrino 
scenarios open a window of opportunity for the search for astrophysical point sources 
of post-GZK UHECR's located at cosmological distances' 3 . In this context, it is 
inter esting to note that AGASA observed a clustering of events on small angular 
(cf. Fig . EJ) - whose statistical significance of occurring higher than 
chance coincidence^^ is still being debatecP™EEBESI however. Intriguingly, the 
integrated flux of cosmic rays contributing to the AGASA event clusters, as shown 
by open squares in Fig. 1101 has a spectrum which is strikingly similar to the one 
expected from cosmogenic neutrinos in a strongly interacting neutrino scenario (cf. 
Fig. (dashed)). Possible correlations of the arrival distributions of UHECR's 
with definite distant astrophysical sources such as compact radio quasard^l in 
particular BL Lacertae objecti^, or GRB's and magnetarJ^2 may give further 
circumstantial evidence for an UHECi' component in EAS data. High statistics data 
from forthcoming cosmic ray facilities such as AugeJ22 and EUS(j25l are requirecP^ 
for these investigations. 



5.1.2. Searches at neutrino telescopes 

The characteristic zenith angle distribution of showers in strongly interacting neu- 
trino scenarios can of course be searched for also at neutrino telescopes where the 
absorbing mat erial is, in addition to the atmosphere, the Earth as well as antartic 
ice (for IceCube 23 ) or water (for ANTARES^J). This is illustrated in Fig. HU (left) , 
which displays the expected zenith angle dis tribution of neutrino-initiated showers 
above 1 EeV in a kilometer-scale detectoJ"3l For Standard Model interactions, the 
distribution (solid curve) is nearly flat for down-going events, and essentially no up- 
going events occur due to very efficient neutrino absorption by the Earth at these 
energies. For models with larger cross-sections, vertical down-going events become 
more frequent, producing more events near cos# zcn ith ~ 1- At zenith angles near 
the horizon, cos# zon ith ~ 0, more of the neutrinos are absorbed and the rate can be 
suppressed. 

Another distinctive observable at neutrino telescopes is the energy spectrum 
of dow n-go ing shower events, which is shown in Fig. ^] (right) for the IceCube 
detectoi^l. The structure of this spectrum can be easily understood. As the 
neutrino energy exceeds the assumed threshold energy of the new interaction (cf. 
Fig-EJj f ne number of events increases dramatically above the Standard model pre- 
diction. Even farther above this energy, however, more of the neutrinos are absorbed 
in the ice before reaching the detector and the event rate is suppressed. This drastic 
"bump" structure in the spectrum indicates the sharply enhanced cross-section at 
the threshold. The peak of this bump occurs at the associated neutrino threshold 



b This window of opportunity is shared with any scenario which exploits neutrinos as "messenger" 
particles, in particular also the Z-burst scenario. 
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Figure 11. Characteristic observables to be studied, in t he context of str ong ly interacting neu- 
trino scenarios, at a kilom eter-scale neutrino detectoJ^I such as IceCubJISI The cosmogenic 
neutrino flux from Refill was as sumed . The dotted lines exploit the cross-section estimate from 
electroweak instantons from Refsl ^lSOl ( cf. Fig. (solid)), whereas the dashed lines represent the 
corresponding estimate from Refs l^ l ^7 l (cf. Fig. [5] (dashed)). The solid lines are the Standard 
Model neutral+charged current predictions. 

Left: Zenith angle distribution of showers generated in neutrino- nucleon interactions. A 1 EeV 
energy threshold for the observed showers has been imposed. 

Right: Energy distribution of down-going showers generated in neutrino-nucleon interactions. 



energy and is mainly generated by charged current electron neutrino interactions. 
The "shoulder" slightly to the left of the bump is from neutral and charged current 
interactions which generate showers less energetic than the incident neutrino. 

Let us ment ion also the possibility to look for enhanced rates for throughgoing 
muons (see RefP^ for a related study) or even spatially compact muon bundled^. 
These signatures, however, rely on details of the final state from the new interaction 
and are, therefore, more model-dependent than the ones discussed above, which 
exploit just generic shower properties. 



5.2. Laboratory tests 

5.2.1. (Quasi-) elastic neutrino- (electron-) nucleon scattering 

As a consequence of dispersion relations, the hypothesized rapid rise in the neutrino- 
nucleon cross-section at large energies is felt in the elastic neutrino-nucleon scat- 
tering amplitude at much lower energieJ^t. Exploiting unitarity and analyticity, 
one may relate the invariant elastic vN amplitudes A± (E) , labeled by the nucleon 



c This has been also pointed out in Refill however in the context of perturbative and model- 
dependent considerations. 
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helicity, with the total vN cross-section via the dispersion relation 

Re, ±(E >- Rc , ± (0> = § P f d*< + «|) , (6, 

where P denotes the principle value of the integral. Suppose the new physics domi- 
nates the neutrino-nucleon dispersion integral © for E' > Eth as hypothesized in a 
strongly interacting neutrino scenario. Assuming that cr"™ is independent of helic- 
ity and energy, and obeys the Pomeranchuk relation cr*^ (E, ±) — crt"^ (E, ±) ^1^5° 0, 
a new contribution to the real part of the amplitude at energy E emerges from (0 , 

Re A ± (E) ~ Re^4±(0) +i- , (7) 

G F /2v / 2 

resulting in 

ReA(£) ncw / ff/100 GeV \ / ggff 
Re^l(£;) S M ~ V^th/10 18 evj ^lOOmb 

for the ratio of the new amplitude to the (perturbative) Standard Model (SM) 
amplitude. Order 100 % effects in the real elastic amplitudes begin to appear already 
at energies seven orders of magnitude below the full realization of the strong cross- 
section. 

Neutrino-nucleon scattering experiments in the laboratory have therefore the op- 
portunity to test strongly interacting neutrino scenarios by searching for enhance- 
ments in the elastic cross-sections. Current experiments at CERN and Fermilab 
reach energies around 100 GeV and, therefore, already start to constrain possible 
scenarios, e.g. p-brane production 1 "- 1 . Elastic and quasi-elastic scattering processes 
may be studied with the help of the HlCSIand ZEU9°-0I detectors at the HERA e ± p 
collider. Its e energy, in the proton's restframe, is around 10 5 GeV and, therefore, 
extends much beyond the energy reach of the above mentioned vN fixed-target 
experiments. However, one-photon exchange dominates the low-energy elastic am- 
plitude for e ± p — ► e ± p to such an extent that the anomalous, new contribution is 
suppressed by a factor of ~ 1/100 compared to (JSJ. On the other hand, possible 
enhancements in the quasi-elastic channels e + p — > v e n and e~p —fD e n, which do 
not suffer from QED dominance, cannot be deduced from model-independent dis- 
persion relations. A separate calculation could be made, however, if certain aspects 
of the new high-energy strong-interaction are assumed. 




5.2.2. Instanton searches 

There is a close analogy^] between electroweak and hard QCD instanton-induced 
processes in deep-inelastic scattering^!. An observably large cross-section for the 
latter processes at HERA is indeed necessary, but not sufficient for an observably 
large cross-section for the former. It seems, moreover, that a > mb electroweak in- 
stanton cross-section necessarily requires that the bulk of inelastic small-Bjorken-a; 
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processes is induced by soft QCD instantons, as has been proposed indeed in RefP. 
Present upper limits on hard QCD instantons from the hJSH and ZEU^SHI 

experi- 
ments are still above the theoretical predictions, but may be improved considerabl 
at HERA II, within this decade. The possible direct production and observatio: 
of electroweak instanton induce d p rocesses in the laboratory will have to wait for 
the commissioning of the VLH 



:cd p 
C 43 . 



6. Summary and conclusions 

We have shown that a simple scenario with a single power-law injection spectrum 
of protons can describe all the features of the UHECR spectrum in the energy 
range 10 17 ^ 21 eV, provided the neutrino-nucleon cross-sections becomes of hadronic 
size at energies above « 10 19 eV. In such a strongly interacting neutrino scenario, 
the cosmogenic neutrinos, which have been produced during proton propagation 
through the CMB, initiate air showers high up in the atmosphere and give thus rise 
to a second, neutrino- induced EAS component, extending well beyond the GZK 
energy. As examples giving rise to the necessary enhancement in u^n, we discussed 
p-brane production in TeV-scale gravity scenarios and Standard Model electroweak 
instanton-induced processes. The model for the proton injection spectrum has few 
parameters from which only two - the power index a and the redshift evolution index 
n - has a strong effect on the final shape of the spectrum. We found that, for certain 
values of a and n, strongly interacting neutrino scenarios are compatible with the 
available observational data from the AGASA and HiRes experiments (combined 
with their predecessor experiments, Fly's Eye and Akeno, respectively) on the 2- 
sigma level (also 1-sigma for HiRes). There are a number of astrophysical source 
candidates, notably neutron stars and GRB's, which may provide the necessary 
conditions to accelerate protons to the required energies, E max > 3 x 10 21 eV, by 
conventional shock acceleration. 

The predicted ultrahigh energy cosmic neutrino component can be experimen- 
tally tested by studying the zenith angle dependence of the events in the range 
10 1S + 20 eV and possible correlations with distant astrophysical sources at cosmic ray 
facilities such as the Pierre Auger Observatory and EUSO, and by looking for bumps 
in neutrino-initiated shower spectra at neutrino telescopes such as ANTARES and 
IceCube. As laboratory tests, one may search for a enhancements in (quasi-)elastic 
lepton-nucleon scattering or for signatures of QCD instanton-induced processes in 
deep-inelastic scattering, e.g. at HERA. 

In summary, strongly interacting neutrino scenarios provide a viable and at- 
tractive solution to the ultrahigh energy cosmic ray puzzle and may be subject to 
various crucial tests in the foreseeable future. 
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